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SUMMARY 

The internal flows of aerospace propulsion en- 
gines have certain common features that are amen- 
able to analysis through computational fluid dy- 
namics (CFD) computer codes. Although the applica- 
tion of CFD to engineering problems in engines has 
been delayed by the complexities associated with 
internal flows, many codes with different capabil- 
ities are now being used as routine design tools. 
This is illustrated by examples taken from the 
aircraft gas turbine engine of flows calculated 
with potential flow, Euler flow, parabolized 
Navier- Stokes, and Navier-Stokes codes. Likely 
future directions of CFD applied to engine flows 
are described, and current barriers to continued 
progress are highlighted. The potential importance 
of the Numerical Aerodynamic Simulator (NAS) to 
resolution of these difficulties is suggested. 


INTRODUCTION 

If an aerospace propulsion system is taken as 
any continuous-flow engine powering a vehicle for 
flight through the atmosphere, then included by 
this definition are gas turbines, ramjets, scram- 
jets and 1 iquid-f uel ea rocket motors. Although 
these engines look and operate very differently, 
they have certain common features that are amen- 
able to analysis, e.g., flow turning, fuel injec- 
tion, mixing, combustion, radiation, etc. In addi- 
tion, certain of them have common flow systems 
(such as labyrinth seals, rotating cavities, cas- 
cade flow, stator -rotor blade interactions, dif- 
fusers, nozzles, etc.) that also can be addressed. 
Propulsion systems, with the exception of rocket 
motors, are usually enclosed in some kind of aero- 
dynamic fairing. Flows over, and into and out of, 
these fairings form a critical part of the whole, 
and provide a link between the external aerody- 
namics of the vehicle and the internal fluid dy- 
namics of its powerplant. It is flows such as 
these with which internal computational fluid dy- 
namics (CFD) is concerned. 

Computational fluid dynamics is emerging as a 
vitally important tool in the design and develop- 
ment of aerospace propulsion systems. It is a tool 
that is being used extensively at present, and its 
continued growth for these applications is assured 
since it makes available to the designer and de- 
velopment engineer information that can be provid- 
ed in no other way. This use and potential of CFU 
is demonstrated in the next section through ex- 
amples taken for convenience, from the gas turbine 
engine. It will be shown that efforts to obtain 
realistic and quantitatively accurate simulations 
are confronted primarily by two technical bar- 
riers. They are: (a) the availability and afforda- 
bility of the computer capacity, and (b) the ac- 


curacy of the physical models that need to be in- 
corporated into the CFD codes. It is anticipated 
that the existence and operation of the Numerical 
Aerodynamic Simulator (NAS) will help to push back 
the first of these barriers and to permit the ex- 
tensive code validation against benchmark experi- 
mental data that eventually will generate improved 
physical models and so, will push back the second 
barrier. 

The application of CFD techniques in the aero- 
space industry has been very successful in exter- 
nal aerodynamics. However, although the status of 
CFD code development is equivalent, successful ap- 
plication to complex internal flows as an engi- 
neering tool has not yet advanced to the same de- 
gree, and the engine manuf acturers appear to lag 
the airframe manufacturers in exploiting CFD by a 
period of at least ten years. It is only in the 
last five years or so when any significant use of 
CFD in powerplant design and development has taken 
place. However, recent progress in the area has 
been rapid. 

Much of the reason for the delay in the use of 
CFD codes as an engineering tool in the design and 
development of propulsion systems can be attribut- 
ed to the considerably more difficult task asso- 
ciated with calculating internal fluid flows. Many 
of the physical simplifications and computational 
shortcuts applicable to problems in external aero- 
dynamics are just not suitable for analyzing in- 
ternal aerodynamic problems. Two-dimensional po- 
tential, and two and three-dimensional Euler equa- 
tion CFD codes, which form the cornerstone ap- 
proach for analyzing external aerodynamic flows, 
are used extensively to analyze turbomachinery 
flows; but in this internal application, they need 
to be calibrated to account for vi scosi ty-generat- 
ed blockage and airfoil trailing edge conditions. 
Boundary layer codes and inviscid/viscous inter- 
action codes are used in conjunction with the po- 
tential and Euler equation solvers to determine 
estimates of airfoil heat transfer and to define 
regions of possible separation. However, at off- 
design conditions where the inviscid/viscous in- 
teractions cause separated regions, this approach 
is not adequate. Parabolized Navier-Stokes codes 
can be used for a class of problems provided that 
regions of separated flow are very small, and 
there exists no significant recirculation of flow. 
There are many situations, however, where the in- 
ternal aerodynamics is dominated by secondary flow 
development, and such large regions of separated 
flow can only be addressed through solution of the 
full Navier-Stokes equations with the attendant 
costs and computer demands. In almost all cases, 
the enclosing geometries are extremely complex. 

This is very much the case for flows in com- 
bustors and dump diffusers, regions of flow on the 
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endwalls near leading edges of compressor and tur- 
bine airfoils, and engine secondary flow systems 
such as rotating cavities and labyrinth seals. 
These problems can only be addressed through the 
solution of the full Navier-Stokes equations. Be- 
cause the bounding geometries in internal aerody- 
namic flows are invariably complex and the flow 
fields are also, large numbers of grid points are 
needed for adequate resolution of gradients; this 
gives rise to high computing demands and the at- 
tendant costs. The perceived difficulties of cal- 
culating internal flows in a realistic and afford- 
able manner did inhibit serious acceptance of CFD 
until such time as computer developments made the 
difficult task feasible. In practice, pragmatic 
engineers in the propulsion field utilize whatever 
solution-approach is physically adequate for their 
purposes and computationally efficient. 

The application of CFD codes in the propulsion 
industry will be demonstrated for a variety of de- 
sign problems associated with the aircraft gas 
turbine engine. It will be shown that some cases 
can be adequately analyzed with the relatively 
simple CFD codes, while other cases require the 
application of a full Navier-Stokes code. Similar 
examples could also be presented for liquid rocket 
motors and ramjet engines. 

The authors wish to thank the Pratt & Whitney 
Engineering Division of United Technologies 
Corporation for permission to present this paper. 
The views expressed are those of the authors. 


CURRENT APPLICATION OF CFD CODES 

In order to illustrate the current state of 
the art, various applications are presented under 
the headings of specific CFD code types. 

Potential Flow Codes 

Since the mid-1970s, subsonic and transonic 
potenti al -f low codes have been widely used at 
Pratt & Whitney to analyze two-dimensional com- 
pressor or turbine airfoil cascades, and axisym- 
metric inlet flow fields. In the last five years, 
there has been an increasing shift toward three- 
dimensional flow analyses. 

One area in which three-dimensional potential 
methods are extensively used is nacelle aerodynam- 
ics. Nacelle-related problems, such as the design 
of inlets and nozzles, involve a combination of 
external and internal aerodynamics. Since the ex- 
ternal flow about the nacelle is essentially irro- 
tational, the potenti al -flow approximation is val- 
id for many of these problems. It is particularly 
applicable to the aerodynamic analysis of inlets 
and fan cowls on turbofan engine installations for 
subsonic transport aircraft. (It is not a very 
good approximation of the physics of nozzle flow 
fields, however, because of the rotational engine 
exhaust flow, strong shocks, and turbulent shear 
layers typical of nozzles at high pressure ratio.) 

Panel methods have seen extensive use in the 
airframe industry for a number of years. They have 
been used at Pratt & Whitney to obtain potential 
solutions for the flow about complex geometries 
such as the installed turbofan nacelle (Lord and 
Zysman, 1986). Most recently, they have been ap- 
plied to the aerodynamic design and analysis of 


Prop-Fan nacelles. The Pratt & Whitney/Allison Gas 
Turbine/Hamilton Standard Prop-Fan propulsion sys- 
tem is a geared counter-rotating pusher. The na- 
celle (figure 1) has an auxiliary air intake for 
the gearbox air/oil heat exchanger. The engine ex- 
haust gas exits the nacelle through discrete lobed 
nozzles located upstream of the propellers. Panel 
codes have proved to be quite useful in analyzing 
the flow into the auxiliary air intake and around 
the lobed nozzles (figure 2). Pratt & Whitney is 
currently using primarily the VSAERO code (Clark, 
et al , 1984). VSAERO is a second-generation panel 
code that exhibits very little leakage, and is 
thus well suited to problems involving combined 
external /i nternal aerodynamics . 



Fan 


LOBED NOZZLE WITH FAIRING 
MACH 0.76 



Figure 2 Distributions of Mach Numbers Over Lobed 
Nozzle of Pusher Prop-Fan - Panel Method 
(Original in Color) * 


The significant advantage of panel codes, of 
course, is that they do not require an off-body 
computational grid. They can handle completely ar- 
bitrary geometries, as opposed to finite-differ- 
ence or finite-volume solvers for which the grids 
and boundary conditions are typically set up for a 
specific component. The computer storage and pro- 
cessor time requirements for panel codes are rela- 
tively modest, well within the capacity of current 
machines. One major drawback is that current 
panel-method solutions are strictly valid only in 
the subsonic, and not in the transonic, flow re- 
g i me . 
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Euler Equation Flow Codes 

Three-dimensional potential methods are not 
generally applicable to propulsion system internal 
flows because of rotational effects (e.g., span- 
wise total temperature and total pressure pro- 
files). This fact helps to explain the intensive 
development effort at Pratt & Whitney of three- 
dimensional (3-D) Euler methods for fan/compressor 
and turbine aerodynamics. 

Euler equation codes are being applied to the 
design and analysis of all the major components of 
the gas turbine engine with the exception of the 
combustor. The particular version used at Pratt & 
Whitney is that developed by Ni (1982). With good 
input/output subroutines, and the incorporation of 
multi-grid methods to speed-up the calculation, 
the Ni Euler-equation solver is extensively used 
by the engineer as just another design tool. The 
Euler code is principally used to determine the 
surface loading distributions on airfoils, end- 
walls and nacelles. It is also used, however, in 

conjunction with 2-D and 3-D boundary layer an- 
alyses to identify potential separated flow re- 
gions and high heat load regions in turbine air- 

foil passages. 

The code was initially developed in two-dimen- 
sional form to calculate the bl ade-to-bl ade pres- 
sure distributions in turbine and compressor air- 
foils. Its main advantage over potential flow 
codes is that it can be applied well into the 

transonic flow regime. The code was extended to 

three-dimensional steady flow for application to 
turbine airfoil passage flows. In addition, be- 
cause of the relatively small viscous/inviscid in- 
teractions that occur in a wel 1 -designed turbine 
airfoil passage away from endwalls, the code has 
proved to be an excellent tool for predicting air- 
foil loadings as a function of span for a variety 
of three-dimensional turbine designs. 

The Ni code demonstrated early that problems 
in existing highly loaded airfoil designs could be 
avoided by changing the airfoil section stacking 
and also the sections to obtain improved pressure 
distributions. In this way, local diffusing re- 
gions could be avoided, and as a consequence, im- 
proved performance invariably was achieved. Some 
examples of this are presented in an American 
Institute of Aeronautics and Astronautics paper by 
Huber, et al (1985). 

One of the biggest problems in applying a 
three-dimensional Euler solver to the design of 
turbomachinery airfoils is the definition of the 
upstream and downstream boundary conditions for 
each airfoil row. An approach adopted to overcome 
this problem was to develop a multi-stage calcula- 
tion procedure, where the upstream boundary condi- 
tion is defined ahead of the first vane and the 
downstream boundary condition is defined behind 
the last blade row. Interface planes between ad- 
jacent blade rows are defined, and flow conditions 
at these planes are calculated as part of the so- 
lution. Implementation of this approach, however, 
requires that certain simplifying assumptions be 
made. Since the flow in actual turbines is highly 
unsteady, and because invariably, for dynamic 
structural reasons, the number of airfoils in 
neighboring rows is different, the flow field at 
the interface planes has to be circumferentially 


averaged before transmitting the information to 
the calculation of the neighboring airfoil row 
which, of course, is in a different frame of ref- 
erence. In principle, the calculation procedure 
could be set up to include an integral number of 
neighboring blades and vanes; and the Euler code, 
being a time-marching scheme, could handle the un- 
steady boundary conditions at the interface plane 
due to the need to change from stationary to mov- 
ing frames of reference between vanes and blades 
and vice versa. However, the current state of the 
art of computers does not allow mesh densities 
that would be sufficiently large to accurately 
capture the flow. It should be emphasized, how- 
ever, that the relatively high levels of turbine 
performance have been obtained with time-averaged 
data and based on the demonstrated accuracy of 
predicting "mean" pressure distributions and flow 
angles. As such, the multi-stage analysis with 
circumferentially averaged interface plane condi- 
tions has provided a significant step in simulat- 
ing the mean flow effects for three-dimensional, 
rotational flow fields in turbines. 

Figure 3 shows a view of a typical two-stage 
high pressure turbine computational mesh. There 
are a total of four airfoil rows contained within 
the computational domain (i.e., two vanes and two 
blades). The lines formed by adjacent airfoil row 
mesh planes are the "interface" boundaries that 
subdivide the computational domain into four sec- 
tors. The three-dimensional flow solution for the 
resultant multisector problem is determined by ap- 
plying the code to all the sectors and imposing an 
interior flow boundary condition to the interface 
mesh points at each time iteration step. The far- 
field boundary conditions for this multistage com- 
putational procedure then become the first vane 
upstream conditions and the second blade down- 
stream conditions. 



Figure 3 Computational Mesh for Turbine Vane/ 
Blade Rows 


Several calculations of two-stage turbines 
have demonstrated the accuracy of predicting air- 
foil and endwall pressure distributions and exit 
flow angles from each airfoil row. In order to 
achieve this accuracy, the secondary flow within 
the passages needs to be reasonably well repre- 
sented. Figure 4 illustrates the predicted vortex 
structure that can exist in the first blade de- 
pending on the flow field exiting the first vane. 
Such secondary flow can have a large effect on the 
flow angles exiting the blade row. 
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Figure 4 Predicted Secondary Flows in First Blade 
Passage - Euler Code 


As the capacity of computers increases, multi- 
stage Euler codes will be extended to simulate in- 
viscid unsteady vane/blade interaction flows, and 
viscous terms will be added to the governing equa- 
tions so that first order viscous effects will be 
simulated in the calculations. 

The Ni Euler code has also been used to com- 
pute the inviscid flow about inlets at angle of 
attack. In this case, a uniform static pressure 
boundary condition is applied at an internal sta- 
tion at or just downstream of the actual fan-face 
station. At low inlet air flows, the flow is isen- 
tropic and the relationship between the captured 
mass flow and the fan-face static pressure may be 
calculated directly from the one-dimensional isen- 
tropic relations. At flow approaching the maximum 
inlet air flow capacity, however, strong internal 
shocks can occur, particularly at high angle of 
attack (figure 5). In the case of nonisentropic 
flow, a static pressure boundary condition is 
specified; but the corresponding actual inlet mass 
flow must be obtained from the computed solution 
by integration of mass flux across an internal 
control surface. A mass-averaged inlet total pres- 
sure loss (shock loss) is also obtained by inte- 
gration of entropy flux across the control sur- 
face. It should be noted that accurate prediction 
of shock loss requires an Euler algorithm in which 
artificial dissipation effects, which lead to 
spurious numerical total pressure loss, are mini- 
mized. 



Figure 5 Mach Number Contours for Engine Inlet 
Showing Shock Formation - Euler Code 


Parameters of interest from a design stand- 
point at high inlet air flow are inlet total pres- 
sure recovery and maximum inlet flow capacity. In- 
let recoveries can be crudely estimated by simply 
adding the shock loss from the Euler calculation 
to the boundary-1 ayer loss obtained from a sepa- 
rate strip boundary-1 ayer calculation. A comparis- 
on with model test data (figure 6) indicates that 
recovery levels are fairly well predicted and that 
the maximum air flow capacity, due to choking at 
the inlet throat, is predicted quite accurately. 



CORRECTED INLET AIR FLOW W c 

Figure 6 Comparison of Predicted and Measured 
Inlet Pressure Recoveries - Euler Code 


Parabolized Navier-Stokes (PNS) Codes 

The Parabolized Navier-Stokes (PNS) approach 
has been applied to both subsonic and supersonic/ 
hypersonic flow fields. At Pratt & Whitney, the 
emphasis has been on development of a PNS code for 
subsonic internal flows in three-dimensional ducts 
of arbitrary cross section and curved centerline. 
A development version of this code is currently 
undergoing extensive calibration. The expected ap- 
plications include fan ducts in turbofan nacelles, 
the radial ly-off set exhaust duct for the lobed 
nozzle in the Prop-Fan nacelle, the circular-to- 
rectangular transition duct in military two-dimen- 
sional nozzle exhaust systems, and several compo- 
nents in the Space Shuttle Main Engine. 

In the PNS formulation, upwind differencing is 
used for transport variables, streamwise diffusion 
terms are neglected, and the streamwise pressure 
gradient is prespecified. In the case of the Pratt 
& Whitney PNS code, global elliptic effects are 
incorporated through use of a coarse-grid Euler 
solution; i.e., it is this solution that provides 
the background pressure gradients. The reduced set 
of PNS equations allows streamwise marching of 
planes of data in space, as opposed to time march- 
ing of three-dimensional data blocks for a full 
Navier-Stokes solver. This results in an order-of- 
magnitude reduction in both computer core memory 
and central processing unit (CPU) time require- 
ments. 

The PNS code has been applied to a well-known 
example case, the Boeing 727 center-engine S-duct 
inlet (Kunz and Rhie 1986). The computational grid 
and a comparison of the Euler and PNS pressure 
distribution with experimental results from model 
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test (Ting, et al, 1975) are shown on figure 7. 
The pressure distributions are here expressed in 
terms of an equivalent Mach number based on up- 
stream total pressure. The inviscid Euler solution 
agrees well with the data over the forward portion 
of the duct but underpredicts the Mach number lev- 
els in the aft portion of the duct. This is as ex- 
pected; the boundary layer builds up with distance 
down the duct and the displacement effect of the 
boundary layer raises the local Mach number lev- 
els. The PNS code viscous solution reflects the 
boundary-layer displacement effect and gives bet- 
ter agreement with the data in the aft portion of 
the duct. 
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SECONDARY FLOW 
AT EXIT PLANE 


Figure 7 Calculation of B 72 7 Center-Engine S-Duct 
Inlet Showing Secondary Flow Development 
- PNS Code 

The predicted secondary flow pattern at the 
exit plane is also illustrated in figure 7. The 
cross-channel pressure gradient induced by the 
second bend drives a circumf erenti al flow in the 
boundary layer toward the top of the duct. A much 
smaller region of secondary circulation is evident 
at the bottom of the exit plane cross section; 
this is the remnant of the secondary flow induced 
by the first bend in the duct. It shows up on a 
total pressure contour as a region of high loss. 

Navier-Stokes Codes 

Flows with dominating recirculation demand 


solution of the full Navier-Stokes equations. Ex- 
amples of such flows are rotating cavities, laby- 
rinth seals and combustion chambers. Two examples 
of such flows are presented. 

Most practical internal fluid flows of inter- 
est are turbulent. This is usually dealt with 
through a simple turbulence management strategy in 
which the Navier-Stokes equations are written in 

Reynolds-averaged or Favre- averaged (density- 
weighted) form and are solved as a stationary 
flow. The turbulent fluxes introduced into the 
equations by the averaging are represented through 
the eddy viscosity hypothesis and gradient diffu- 
sion. Turbulence modeling is used to find rela- 
tionships for the eddy viscosity. This modeling 
may be appropriate for both high and low turbu- 
lence Reynolds numbers, or a high Reynolds number 

model may be used together with turbulent wall 

functions for the near-wall region. For variable 
density and chemically reacting flows, the energy 
equation and species transport equations are in- 
cluded. The partial differential equations are ar- 
ranged into a general form that consists of con- 
vection and diffusion terms, and source terms de- 
scribing the generation and dissipation of the de- 
pendent variable. A heat release rate expression 
is provided for reacting flows. Finite differenc- 
ing is used to discretize the equations. Since 
pressure does not appear directly in the momentum 
equations or the continuity equations but is a de- 
pendent variable nonetheless, an algorithm is used 
to obtain a pressure field and correct the veloci- 
ty field in such a manner that the continuity and 
momentum equations are simultaneously satisfied. 
An implicit solution algorithm is used, with the 
initial guesses for the field variables being 
iteratively updated until convergence is reached 
(Sturgess, 1983). The liquid fuel is accounted for 
through a Lagrangian spray model that is coupled 
with the Eulerian gas flow field through a par- 
ticle-source-in-cell technique (Sturgess, et al , 
1985b). 

Part of the difficulty of presenting internal 
flows is that they are totally enclosed. In addi- 
tion, the enclosing geometry can be extremely com- 
plex. The combustion chamber of the gas turbine 
engine is one such component; this is illustrated 
on figure 8 which is a view ot an annular combus- 
tor presented from an interior perspective. Al- 
though the flow can be considered as periodic 
about fuel injectors, it may be appreciated that 
definition of the many streams (film cooling air, 
combustion air jets, dilution air jets, dome cool- 
ing air, swirler air, and liquid fuel), entering 
the computational domain, as well as the basic 
bounding geometrical contours, demand a large com- 
puter resource. Even when reduced to a repeating 
segment, the combustor geometry presently has to 
be greatly simplified, an example of which is 
given in figure 9. 

Calculation of the simplified combustor of 
figure 9 yields a tremendous amount of useful in- 
formation concerning the details of the reacting 
flow contained therein. Figure 10 gives an iso- 
therm contour plot at a cross-section in the dome 
and shows details of the film cooling, the dome 
cooling associated with the heatshield, and the 
results of local air admission around the fuel in- 
jector. Note that the temperature field is not 
symmetric about the axisymmetric fuel injector and 
heatshield geometries; this is due to the swirling 
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air introduced by the airblast fuel injector. An 
additional cross-section is contained in figure 
11. This section is taken at a plane through the 
second row of air-addition ports. Streaklines are 
provided as motion-cues and are colored with the 
appropriate color for the local temperature field. 
The complex character of the flow field can be 
appreciated. 



Figure 8 Interior View of an Annular Combustor 
for a Gas Turbine Engine 
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Figure 10 Isotherm Contour Plot Close to Dome in 
Simplified Dome - NS Code (Original in 
Color) 
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Figure 11 Combined Streakl ine/Temperature Plot at 
Cross-Section Through Second Row of 
Air-Ports (Original, in Color) 

Figure 11 demonstrates the dominant character- 
istic of combustor flows, i.e., they are wholly 
elliptic. This is best illustrated by means of an- 
imated streaklines, of which figure 12 is a 
"still" rendering. 



Figure 12 Isometric Projection of "Frozen" 
Animated Streaklines in a Model Primary 
Zone (Original in Color) 

Calculations in the combustor such as those 
shown are most useful to understand the flow pro- 
cesses, for determining the origin of hot-spots on 
the confining liners and dome (Sturgess, 1980), 
and in investigating geometric changes to develop 
the outlet temperature distributions. This is all 
done by directly viewing the reacting, hot-gas 
flow field, rather than by measuring its effects 
on the confining boundaries and then postulating a 
field which is to be modified in some way. 

The calculations shown were performed with a 
Cartesian coordinate system on a 50 x 40 x 41 
grid, requiring a storage of 30 megabytes and 
about 11 CPU hours on an IBM 3090 (fast scalar) 
computer to reach a somewhat arbitrary convergence 
level of 5 percent (residual source sum). At this 
level, the definition of the combustor is inade- 
quate, and calculations at the 200 megabyte level 
are being explored, although it is estimated that 
300 megabytes might be necessary for adequate 
representat ion of the geometry and resolution of 
the flow gradients. The solutions currently ob- 
tained cannot be considered to be grid- independ- 
ent. Thus, when calculations like this are used 
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for parametric diagnostic studies, all variations 
contemplated must be calculated on the same grid; 
this necessitates very careful thought and ad- 
vanced planning in establishing the initial grid. 

The accuracy of combustor simulations is being 
determined at present not by the numerics or the 
physical modeling but by computer capacity and the 
cost of solution. 

Prepared originally for the combustion cham- 
ber, the availability of a general Navier-Stokes 
code written in modular form containing a suite of 
switchable physical models and an arbitrary geome- 
try capability, enables application to be made to 
internal flows with recirculation that arise in 
many areas of the gas turbine engine, including 
secondary systems as well as the main gas path. An 
example of such an important secondary system is 
the internal cooling of turbine aerofoils where 
there is a need to understand and predict the ef- 
fects of rotation on local heat transfer and pres- 
sure loss, especially for blades that have multi- 
pass coolant passages. Accurate calculation of 
these quantities is very important as it can lead 
to an improvement in the reliability of predicting 
turbine airfoil temperatures and, ultimately, 
blade life. Accurate estimates of temperature and 
life permit available cooling air to be used most 
effectively and reduces its negative impact on 
turbine performance. 

In order to demonstrate the potential for 
this, a simplified geometry of a blade passage was 
calculated, as shown on figure 13 (Sturgess and 
Oatta, 1987). It consisted of a single outflow ley 
and an inflow leg. The grid was a coarse 33 x 11 x 
5, and these demonstration calculations were not 
considered to be grid independent, and the resolu- 
tion is low. For the baseline calculation of 600 
RPM, the Reynolds number was 30,000 and the Rossby 
number was 0.174; the passage aspect ratio was 
unity and its length-to-width ratio was 12. 



Figure 13 Schematic of Simulated Turbine Blade 
Coolant Passage 


Figure 14 gives the flow field development 
along the passages at 600 RPM and at the two 


cross-sections shown in figure 13 for several 
rotational speeds. The flow visualization tech- 
nique is by streaklines, and results for four 
rotational speeds of 0, 60, 600 and 1 900 RPM are 
presented. The development of secondary flows due 
to the influence of Coriolis forces can be seen. 
In the outflow (away from the axis of rotation) 
leg of the passage, a pair of counter-rotating 
vortices develop. The vortex centers are shifted 
slightly from the centerlines of the leg towards 
the pressure side, but the size and strengths of 
the vortices are about equal. As the vortex pair 
enters the turn, their rotational velocity is 
overcome, and the flow changes direction as mass 
is forced to the pressure side of the turn. The 
flow entering the passage inflow leg from the turn 
is therefore forced into a right-angled corner. As 
it escapes from the corner to begin flowing down 
the inflow leg (towards the axis of rotation), the 
air has no choice other than to establish a single 
vortex. The action of the turn is thus to coalesce 
the vortex pair formed by Coriolis forces into a 
single vortex completely filling the passage in- 
flow leg. Other sections of the passage show that 
the the double vortex structure of the original 
system begins to reestabl i $h itself as the flow 
proceeds down the inflow leg towards the axis of 
rotation, although initially this vortex pair is 
not symmetrical, of course. The flow behavior ob- 
served is intuitively correct. 
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Figure 14 Flow Visualization in Simulated Coolant 
Passage for a Range of Rotational 
Speeds - NS Code 
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Although the code can give adequate quantita- 
tive accuracy in three-dimensional rotating pas- 
sages of very simple geometry, as was established 
by subsequent calculations (Sturgess and Datta, 
1987), the geometries of practical interest re- 
quire more computer capacity than is currently 
available. For example, the use of trip-strips in 
the passages to enhance heat transfer requires 
that the induced disruption of the laminar sub- 
layer in between strips be both calculated and re- 
solved in detail. Such a calculation demands the 
use of appropriate low Reynolds number turbulence 
models and many grid lines adjacent to the walls. 
Again, progress is limited by the available com- 
puter f aci 1 ity. 


FUTURE DIRECTION 

It has been shown by way of example, how cur- 
rent CFD codes are being used in the design and 
development of aerospace propulsion systems. In 
conjunction with highly refined design criteria, 
advanced materials and advanced manufacturing 
technology, gas turbine component performances 
have progressively increased so that the thrust 
specific fuel consumption (TSFC) of commercial and 
transport gas turbines has been reduced at a rate 
better than 1% per year over the past two decades 
(figure 15). Also shown in figure 15 is that, with 
the introduction of advanced propulsion concepts 
like the Prop-Fan or the ducted fan, the potential 
for continued significant TSFC reductions will 
continue. Improved modeling of the gas path flows 
using the developing CFD codes in the new propul- 
sion systems and in advanced versions of existing 
gas turbine engines will play an important role in 
achieving the full potential of these propulsion 
systems. 
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Figure 15 Improvement in Gas Turbine Engine Fuel 
Consumption with Time 

For military gas turbine engines, high thrust- 
to-weight ratio and specific thrust are major re- 
quirements. Figure 16 shows the improvements in 
thrust-to-weight ratio that have been made in re- 
cent years, and also a projection of what is re- 
quired in the future. It can be seen that the rate 
of improvement over the next 1 5 years must be 
about twice the current rate. It has been identi- 
fied that approximately half of the improvements 
can be achieved with advanced materials, and the 
other half needs to be achieved by improved com- 
ponent efficiencies and by reducing the air "leak- 
age 11 from the propulsion stream. Once again ad- 
vanced CFD codes will be needed to achieve the 
necessary component performance goals and to en- 
sure that as much as is possible of the air enter- 
ing the engine goes to generate useful thrust. 

The successful applications of CFD codes ob- 




JT9D 

IMPROVEMENTS 


FUTURE 

• 3-D VISCOUS AERO o 

► IMPROVED COOLING 

► IMPROVED MATERIALS 

I I 1 


o) 


PW4000 


70 


1980 


1990 


2000 


tained to date should not give rise to compla- 
cency. To the contrary, if the aggressive military 
and commercial goals needed to keep this nation 
competitive with outside threats are to be met, 
then significant advances in the development and 
application of CFD codes will need to be made. Two 
major limitations exist that currently are imped- 
ing further progress in the development and use of 
CFD codes. They are: the inadequacy of physical 
models that need to be incorporated in order to 
simulate realistic engine flow fields, and the 
limitations in computer capability associated with 
speed, memory capacity and cost of calculation. 



Figure 16 Increase in Engine Thrust to Weight 
Ratio with Time 


The first and perhaps the most obvious defi- 
ciency in current CFD codes is the lack of realis- 
tic physical models. Errpir leal data from many dif- 
ferent types of test are used to calibrate the 
codes for application to gas turbine design. For 
example, in turbomachi nery flows, surface shear 
forces and trailing edge conditions are establish- 
ed in order to obtain the required accuracy in 
predicting airfoil loadings, air angles and weight 
flow. Turbulence models that are necessary to pro- 
vide closure to the Reynolds-averaged Navier- 
Stokes equations have evolved to the point where 
relatively simple flows can be simulated, but they 
are insufficient to be used as a reliable means 
for predicting airfoil loss and heat transfer. For 
the reasonably near future, it seems unlikely that 
the complex flows in airfoil boundary layers and 
in free shear flows will be calculable. Thus, if 
CFD codes are to be used in the future to optimize 
on minimum loss or heat transfer, a large input 
from key experiments will be required. Similarly, 
the interactions between turbulence and heat re- 
lease in the combustor are not well understood and 
can only be modeled in crude fashion, for which 
empirical input is required. The problems in gen- 
erating these data are twofold: what are the key 
experiments to be performed that are needed to de- 
velop codes which can be applied to flows in an 
engine environment, and how will the large number 
of these experiments be funded? Attempts at an- 
swering these questions are beyond the scope of 
this paper, but they need to be addressed by all 
who are in the business of developing and applying 
CFD codes. 

No discussion of the limitations of physical 
modeling can omit the fact that almost all of the 
CFD codes used in gas turbine design are steady, 
and that the so-called Navier- Stokes codes usually 
do not solve the Navier-Stokes equations, but em- 
ploy a statistical turbulence management strategy 
and solve the time-averaged Reynolds or Favre 
equations. Yet, of course, the flow in most gas 
turbine components is highly unsteady. In fact, it 
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is quite remarkable that so much progress in cal- 
culating component flows has been made without 
considering unsteadiness as a first order factor. 
Again, using an example from turbomachinery flows, 
the flow in a specific turbine or compressor air- 
foil passage is very much dependent on the flow 
field exiting the upstream airfoil row. As dis- 
cussed in the section on Euler codes, because the 
upstream row is moving relative to the specific 
airfoil passage in question, the flow field from 
the upstream row currently is pitchwise averaged 
and the flow is assumed steady and circumferen- 
tially constant as it enters the airfoil passage 
that is being analyzed. From the work of Sharma, 
et al (1985) it can be shown that this simplifica- 
tion results in the suppression of vorticity en- 
tering the airfoil passage. Since the function of 
the airfoil passage is to turn the flow (and the 
turning angle can be over 100 degrees for tur- 
bines), the redistribution of the entering vorti- 
city is instrumental in determining the exit flow 
angles. Thus, if much of the inf low-vorticity is 
suppressed, it is not possible to obtain reliable 
exit flow predictions from an essentially steady 
CFD code. Similarly, in the combustor, the current 
calculation of the distribution of fuel across the 
reaction zone depends on the trajectories of the 
fuel droplets (which should include the effects of 
turbulent dispersion of the droplets and modula- 
tion of the gas-phase turbulence by the presence 
of the droplets) and the assumption of gradient 
diffusion. However, observations of the flow pat- 
terns in combustor primary zones show that nonsta- 
tionary effects can be dominant. Thus, the actual 
mixing of fuel and air is largely by processes 
that are not accounted for in the calculation pro- 
cedure which solves time-averaged equations. One 
solution to these problems is to develop unsteady 
three-dimensional CFD codes, and starts have been 
made in this direction. Such codes introduce such 
a level of complexity and calculation cost, how- 
ever, that it is unlikely that unsteady codes will 
be used significantly in the design of gas turbine 
components for several years (Sturgess, 1984). 

The other significant limitation in developing 
and applying CFD codes is the current state of the 
art of computer hardware. In gas turbines, the 
geometries and the flow fields are highly complex; 
thus, large numbers of grid points are desired by 
the code developers in order to predict the flow 
field with the accuracy required by the designer. 
The designer also wants to obtain answers fast, 
make some changes and re-run the code until appro- 
priate design criteria are met. Since designs are 
generated within a given budget, the designer also 
wants the cost of each run to be reasonable in or- 
der to stay within that budget. The code developer 
of course also wants inexpensive and fast comput- 
ing capability so that the parametric studies nec- 
essary for calibrating the code for a wide range 
of geometries and flow conditions can be conduct- 
ed. Many papers in the literature have realistic- 
ally assessed the computer hardware needs for the 
future, e.g., Sturgess (1985a) for the combustor, 
and the Numerical Aerodynamic Simulator (NAS) has 
a well-established program in place to meet many 
of these needs. Within this context, emphasis must 
be placed on re-configuring CFD codes to be opti- 
mally efficient on vector ana parallel-processing 
computer archi tectures. Emphasis must also be 

placed on exploiting such techniques as adaptive 
grid procedures to reduce as much as possible, the 
computer resources needed for solution of a spe- 


cific problem. 


CONCLUDING REMARKS 

It has been shown that internal flows for 
aerospace propulsion are now being extensively 
calculated in the industry as part of the design 
and development procedure of engine components. 
Although the illustrative examples presented have 
been taken exclusively from the gas turbine en- 
gine, this is so for ramjets, scramjets and 
1 iqu id-fueled rocket motors also. A variety of 
codes is used, and the physics embodied in these 
codes range from potential flow to the Navier- 
Stokes equations. Selection of the appropriate 
code depends on the component flow required to be 
calculated and on the objectives required of the 
calculation. The codes are currently making a use- 
ful and important contribution to achieving the 
required performances of the various engines. 

While the utility of the current generation of 
codes cannot be denied, there are some limitations 
to the accuracy of the calculations produced, and 
these were briefly touched upon. These limitations 
can be succinctly described as: (1) lack of physi- 
cal realism in the modeling, and (2) restraints 
imposed by available computer resources. 

The Numerical Aerodynamic Simulator (NAS) rep- 
resents a great resource with immense potential 
for assisting with the calculation of internal 
flows in aerospace propulsion engines of all 
types. It is an opportunity for the present bar- 
riers to accuracy in such calculations to be over- 
come. The NAS, because of its high speed and large 
storage capacity, can be used as an experimental 
wind tunnel, or in the present circumstances, as 
an experimental engine component. In this role, it 
can be used to make calculations that are free 
from computer-resource limitations; therefore, it 
can be used to improve the understanding and mod- 
eling of relevant physical phenomena and to pro- 
vide calibration information that is necessary for 
simpler, design-oriented codes intended for regu- 
lar engineering use. 

For example, the tip clearance flow field is a 
major contributor to loss in compressors; but be- 
cause of the prohibitive number of grid points 
which would be necessary to incorporate it accu- 
rately in a numerical treatment of a full blade 
row, it is currently represented by very crude em- 
pirical modeling. The NAS would permit calculation 
of the actual tip clearance flow field with suf- 
ficient resolution when embedded in a full three- 
dimensional representation of blade row geometry 
to lead to more fundamental models. 

In the case of non-stationary flows, the capa- 
bilities of the NAS can permit the use of unsteady 
codes for relevant cyclical periods. For turbo- 
machinery applications, blade and vane row inter- 
actions can be examined and methods developed to 
model average flows at computational interfaces. 
This information can be incorporated in steady 
flow design codes, which will then more accurately 
model the interaction effects. Similarly, non-sta- 
tionary mixing and chemical reaction processes in 
the combustor could be accurately modeled on the 
NAS facility, permitting full evaluation of the 
relative magnitude of the effect of time-dependent 
processes on the time-average solution. 


The examples given above suggest how the NAS 
might contribute to engineering design capability. 
How useful the facility actually is will depend on 
the aerospace propulsion industry learning how to 
exploit its capabilities to the best advantage. 
Certainly, it seems apparent that the computa- 
tional capacity offered by the NAS will be as nec- 
essary to achieving the goals of the new engine 
for the Advanced Tactical Fighter (Petty, et al, 
1986) and the lightweight stoichiometric engine, 
as it is to the Space Shuttle Main Engine and the 
aerospace plane. 
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